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We present the exact Bethe Ansatz solution of a multichannel model of one-dimensional correlated 
electrons coupled antiferromagnetically to a magnetic impurity of arbitrary spin S. The solution 
reveals that interactions in the bulk make the magnetic impurity drive both spin and charge fluc- 
tuations, producing a mixed valence ntmp ^ at the impurity site, with an associated effective 
spin Seff = 5-1- |nimp|/2 in the presence of a magnetic field. The screening of the impurity is 
controlled by the size of the impurity spin independently of the number of channels, in contrast to 
00 ' the multichannel Kondo effect for free electrons. 
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PACS. 75.20.Hr - Local moment in compounds and alloys; Kondo effect, valence fluctuations, 
D ' heavy fermions. 71.27.+a - Strongly correlated electron systems. 

^ ' Introduction. - Impurities play an essential role in correlated electron systems, in particular in ID where even a 
small amount of defects may drastically change the properties of the system. While the effect of a single potential 
scatterer is by now fairly well understood 0], the case of a dynamical scatterer in a ID interacting electron system - 
^ \ like a magnetic impurity - largely remains an open problem |^ . 

r ' Of particular interest is to understand what happens when the electrons are coupled to a local magnetic moment 
^ \ (of magnitude S) in several degenerate channels m. In the case of a 3D Fermi liquid, the impurity induces strong 
Q ■ correlations among the electrons, with the low-temperature physics depending on the relation between the impurity 
O ', spin and the number of electron channels |3|-|^] . It is natural to ask whether the correlations inherent in a ID electron 
system (Luttinger liquid, |^) will modify this behavior. Does novel effects appear, or does one recover the same 
' multichannel Kondo physics as for free electrons in 3D? Apart from the possible experimental relevance of these 
^ ' questions to the study of artificial impurities in mesoscopic devices their resolution is an interesting issue in its 
QO \ own right. A study is also motivated by the recent interest in the Kondo effect in the high- Tc cuprates jl^: it is 
well known that at least two orbitals, 3dx'^^y2 and Sd^^, play an essential role there, so the Kondo effect should have 
multichannel nature. An analysis of the simpler, analog problem in ID - with correlated electrons - may provide 
important clues for how to model the effect in the cuprates. 
t — In this Letter we attack the problem by considering an integrable model of a magnetic impurity embedded into a 

0^ multichannel interacting electron system. We here take a multichannel extension of the supersymmetric t—J model in 
ID with the electrons coupled to a localized spin (of magnitude S) by an antiferromagnetic exchange interaction 
that preserves integrability. This makes applicable a Bethe Ansatz approach, allowing us to obtain exact results 
C ' for the groundstate properties as well as the finite temperature behavior. Our solution reveals that there are two 
. distinct processes that govern the zero temperature response due to the impurity: First, the singlet Cooper-like pairs 
P5 ' of electrons present in the ground state break up and get temporarily trapped when scattering off the impurity, thus 
O producing an impurity valence riimp with an associated effective impurity spin Seff = "S*-!- |nimp|/2 in the presence 
of a magnetic field H. Secondly, this effective composite spin gets screened by unbound electrons excited by H from 
the spin singlet groundstate. The type of low-temperature behavior that emerges depends only on the size of the bare 
impurity spin S, and is insensitive to the number of channels. This property is different from that of the ordinary 
multichannel Kondo effect in a Fermi liquid. 
d [ The model. - The Hamiltonian Hhost of the multichannel extension of the supersymmetric t — J model can be 
written as Hhost = J2j '^j.j+i^ ^'^^^ 

Here j — 1, L is a site index, cr = ± denotes the spin projection, / = 1, . . . , m is a flavor quantum number indexing 
the available electron channels, and P is a projector on the subspace of singly occupied states. All indices are summed 
over. The scattering matrix for two electrons is given by X{pi—pj) = [{pi—pj)I+iP'']^[{pi—pj)I — iP-f/[{pi—pj)'^ + l] 
with pij the corresponding charge rapidities, and P'' and P-^ the permutation operators in spin and flavor subspace, 
respectively. The X-matrices satisfy the Yang-Baxter equation X^^{pi — p2)X^^{pi — P3)X^^{p2 — Ps) = X^^{p2 — 
P3)X^^{pi — p3)X^^{pi —P2), thus ensuring the integrability of the model. 
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We now insert an additional site on the lattice, call it 0, attach a magnetic impurity to it, and couple it to 
the electrons on neighboring sites (sites 1 and L, given periodic boundary conditions). To preserve integrability 
this interaction must be judiciously chosen. We here follow the strategy pioneered in and define the impurity- 
host interaction via an electron- impurity S-matrix S that satisfies the Yang-Baxter equation X^^(pi — p2)S^'^{pi — 
Po)S^'^{p2 —Po) = S'^^(jp2 — pq)S^^{pi — po)X^^{pi — P2), with po measuring the impurity energy level. This approach 
assures that the model remains integrable in presence of the impurity. The S-matrix can still be chosen in a number 
of ways , and here we take it to be similar to that of the multichannel Kondo problem in a free electron gas 0| 
and let it act nontrivially only in the spin subspace. Writing out the components, 

p ~\~ ^ 

with a'' the Pauli matrices {k — x,y,z), S'' the impurity spin matrices, \M\ < S the component of the impurity spin 
S (with unprimcd/primed indices referring to incoming /outgoing states), and a{p) = [(p^ + l)/(p^ + {S + 1/2)^)] 1/2, 
The corresponding impurity Hamiltonian Tiimp can be written in the form: 

Ti-imp = J{HLfi + Hos + {'Wl.o, ^0,1}) + (1 ^ 35(5 + 1) J - —)Hl,i 

+ 2poJ[{nLM + noA),nLA] . (3) 

Here Tio.i and TLl.o a-re generalized permutation operators of a particle with spin 5 but no flavor (impurity) and an 
electron (carrying both spin and flavor), while Hl,i is defined in Eq. (1). The function J = [pq -I- (5 -I- 5)^]"^ plays 
the role of an effective exchange constant between the impurity site and neighboring sites. We should mention that 
the structure of Himp for the periodic chain in (|^) simplifies considerably when the impurity is located at the edge of 
an open chain, with the impurity connected to the host via only one link, with coupling constant J. It is here worth 
pointing out that the low-energy behavior of impurities in integrable open and periodic chains of correlated electrons 
are qualitatively the same, as was shown recently in We should also point out that our approach is different 

from that recently advocated by Wang and Voit in their study of the ferromagnetic Kondo effect in a Luttingcr liquid, 
where the impurity is simulated by a boundary potential [p^ . In contrast, in our formulation the full dynamics of the 
impurity is retained. 

Bethe Ansatz equations. - The model thus constructed can be diagonalized exactly by an algebraic Bethe ansatz 
pGf . The procedure is rather cumbersome, and we here only give the result. The eigenstates are characterized by m-|-2 
sets of quantum numbers: The charge rapidities {pj}^^i (with N the total number of electrons), the spin rapidities 

{^a}a=i (-^ counting the number of spin-down electrons), and the m sets of flavor rapidities {/i/jjjjtli , i = 1, 2, m 
(with m^*^ = X]feLi+i n'^^\ n'^''^ counting the number of electrons in channel A;, fc = 1, 2, m). Each state corresponds 
to a particular solution of the nested Bethe Ansatz equations 

n "ii ' - 4^^') = n - <') ' 

T=± (3=1 7=1 

[e2S+i(Pi -Po)e2(Po -Pj)Y''^e.[{pj) = J| ei{pj - X^) J| eiip-^^'^ - Pj) 

a=l /3=1 

JV J\/ 

e2s(AQ -po) ]^ei(AQ - p^) = e2(Aa - Xp) , (4) 
j=i 0=1 

where e„(a;) = {2x + in) / {2x — in) , pL"^^ — pj, m^^^ — N, m^"^'^^'^ — 0, and L is the number of lattice sites (not counting 
the impurity site). We shall assume that no external fields couple to the electron channels, and can hence confine our 
attention to the flavor-singlet subspace. The energy (up to an additive constant) and the magnetic moment are then 
equal to E ^ (P^j + (1/4))"^ and 5^ = A^/2 + 5 - M, respectively. 

The thermodynamics of the model is described by the usual string hypothesis |p^ . In the thermodynamic limit, 
with L,N,M,m'^^^ — > 00, their ratios being fixed, the model has the following possible excitations: (i) unbound 
electrons with charge rapidities pj] (ii) spin singlet Cooper- like pairs with pj = Xj ± i/2; (iii) spin strings (bound 
states of any number of "down spins"); and (iv) flavor strings for each channel. (Since the host interactions in (|l|) carry 
opposite signs in spin- and flavor subspaces, bound states of different channels as well as bound states of electrons 
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and spin strings are suppressed.) Introducing distribution functions (densities) for particles and holes for each class 
of excitations we can write down the corresponding thermodynamic Bethe ansatz equations. Then, by minimizing 
the free energy, we extract the integral equations for each excitation class. These equations have the same structure 
as those of the multichannel Kondo problem for free electrons JT^ , and differ only in the driving terms (which are 
independent of energy and density). This implies that in the high-temperature limit, where the driving terms are 
unimportant, the effect of our impurity is similar to that of a Kondo impurity in a multichannel free electron gas. 

Groundstate properties. - Let us focus instead on the groundstate properties of the impurity, and study how it 
depends on the host band filling and an applied magnetic field. In the zero temperature limit, T — > 0, only excitations 
of classes (i), (ii), and (iv) (for unbound flavorons and pairs of them, i.e. flavor strings of length 1 and 2 for each 
channel) can have negative energies. Thus, the filling of these states is determined by the Dirac seas of the groundstate 
of the model. In the singlet flavor sector we can solve the equations for densities and energies of the flavorons as 
functions of the densities and energies of unbound electrons and singlet spin-charge pairs. As a result, we obtain the 
ground state equations for unbound electrons and pairs: 

Ph{p) + {l-ai*si)'k{p{p)+ai'ka{\))=ai{p) + ^[a2s+i-a2]ip-Po) , 
cr,i(A) + (1 + 02) * (1 - ai ★si) ★ (cr(A) + s*p(p)) = a2(A)+2x[a2S+2-a2S-a3-ai](A-po) , (5) 

where * denotes convolution, p{ph) and a{(7h) are densities for unbound electrons and pairs, respectively. The 
integration intervals are given by \p\ > B and |A| > Q, with Q and B playing the role of Fermi points for unbound 
electrons and pairs, respectively. The Fourier transforms of the kernels a„, s and si are given by exp(— n|u;|/2), 
cosh^^ /2uj, and sinh[(m — l)|w|/2]/ sinh(TO|ci;|/2), respectively. Eqs. (|^) are linear in the densities, and the driving 
terms of the host and the impurity are additive. Separating the densities into bulk and impurity parts then allows 
us to calculate the valence and the magnetization of the impurity in the groundstate. In the absence of an external 
magnetic field we have B ^ 00 (no unbound electrons). The limit Q 00 corresponds to a vanishing pair density, 
while Q — > is the limit of l/2TO-fillcd bands of the host (no pair holes, corresponding to a vanishing Fermi velocity). 
For the nonmagnetic groundstate it follows that the effective valence of the impurity varies as a function of electron 
number from n.imp — for vanishing pair density to nimp ~ —rii for l/2TO-filled bands, with the valence measured w.r.t. 
the groundstate of spin-paired electrons (i.e. its negative sign indicates an excess of pair holes due to the presence of 
the impurity). This unusual behavior is caused by the host interaction in (Q) which supports a groundstate with a 
Dirac sea of bound singlet Cooper-like electron pairs, not present in the free electron gas. 

Of special interest is the behavior of the impurity magnetization 5*^?^^. With no magnetization in the bulk, we have 
^imp ~ ^- By turning on a magnetic field, the number of unbound electrons increases while the number of singlet 
pairs decreases (as required by electron number conservation). By eliminating the pair density cr(A) from the second 
equation in (|^), we obtain the Fredholm equation 

Phip) + pip)-F-kp = s*(7h + s(j)) + jj-s-ka2sip~Po) , (6) 

where F{uj) = 1 — tanh(|w|/2)[l — exp(— toIwI)]"^, and the integration over the pair hole density is over the finite 
interval [~Q,Q]. This yields an explicit connection between the densities of unbound particles (electrons or holes) 
and spin singlet pair holes. As the Zeeman splitting is typically much smaller than the Fermi energy, we can neglect 
the influence of the pairs on the impurity magnetization as long as the magnetic field is sufficiently weak. For this 
case {H <^ 1), choosing Q = (i.e. l/2TO-filled bands) and assuming that \po\ ^ B, we can solve Eq. (^) exactly by 
reducing it to a sequence of coupled Wiencr-Hopf integral equations. In this way we obtain two distinct regimes for 
the behavior of S^^^^ with magnetic field. If 5 > 1/2 the impurity spin becomes asymptotically free: 

= m[1 ± f (I log(ff/TH)|)-i + ...], 5 > i (7) 

with Th = 27r(|^)'"/2TK/r(f ), and Tk oc exp(-7rpo) playing the role of the Kondo temperature. When H < Th, 
fi — S, and the upper sign in (|^) gets selected. On the other hand, with Th H ^ 1, p, = S + m/2 and the 
lower sign in (|^) is selected. Note that for this case H is still much smaller than the spin saturation field in the bulk 
(corresponding to a transition to a ferromagnetic bulk state). Also note that the underscreened behavior in ^ for 
H « Th appears only for m < 25 in a free electron gas 0, in contrast to the correlated host studied here. For the 
exactly screened case S = 1/2, 
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producing a finite zero-field magnetic susceptibility. Again, this is different from a free electron host where this 
behavior sets in for m = 2S. Most importantly, the case of overscreening (rn > 2S for a free host, with critical scaling 
of the impurity magnetization) does not exist in this interacting electron system. Thus, electron correlations in the 
host due to direct electron-electron interaction here suppresses the critical overscreened behavior of a Kondo-like 
impurity, similar to what happens in the presence of a channel anisotropy for free electrons . It is important to 
realize that by increasing the magnetic field we deplete the number of singlet pairs. Therefore, as we have already seen, 
by varying the field we can smoothly tune the effective impurity valence riimp from — m to 0. In the underscreened case 
this means that the uncompensated part of the effective impurity spin Seff = 5"-!- |«imp|/2 decreases with increasing 
field, and recovers its zero-field value S at the spin-saturation field. 

Low-temperature thermodynamics. - We now turn to the low-temperature thermodynamics. As we have already 
mentioned, the thermodynamic Bethe ansatz equations differ from those for the free electron gas multichannel Kondo 
problem only in the driving terms, and we can hence use an analysis similar to that in 19 1. The main idea is to 



rewrite the T ^ thermodynamic Bethe ansatz equations for all excitations on universal form. This can be realized 
in the "Kondo limit" (with suppressed charge fluctuations), yielding for the universal energy potentials (pj 19|: 



0,(A) = .s*log(l + e^-0(l + e'^^+0-'5,,„e^^ , (9) 

with limj_+oo(</'j/j) = H/T. It follows that in this limit the impurity low-temperature free energy can be written as 
a function of (pm ■ 

p _ ^ r rfAa2s(A)log(l + e^-) 
™^ " y_oo 2 cosh[7rA + \og{T/TK)] ' ^'"^ 

The qualitative behavior of Fimp is independent of the relative values of m and S, in analogy with the groundstate 
properties. In the underscreened case S > 1/2, the impurity magnetic susceptibility Ximp shows Curie-like temperature 
dependence while the specific heat Cimp exhibits a Schottky anomaly at T ^ H and a Kondo resonance. For the 
exactly screened case 5' = 1/2 local Fermi liquid behavior holds: Ximp ~ 'm/2'KTK, Cimp ~ i:mT /(m -\- 2)Tk- 

Discussion. - To conclude, we have found an exact Bethe Ansatz solution of the problem of a Kondo-like antifer- 
romagnetic impurity embedded into the multichannel supersymmetric t — J model. Our solution reveals that there 
are two mechanisms governing the magnetic behavior of the impurity: holes in the distribution of Cooper-like pairs 
increase the effective spin of the impurity, while unbound electrons screen this effective value. Note, that for the t—J 
model there is no spin gap, so the unbound electrons and Cooper-like pairs coexist even for small external magnetic 
fields. In contrast to the multichannel Kondo effect in a Fermi liquid, the screening with unbound electrons are 
realized in only two different ways independent of the number of channels, and controlled only by the value of the bare 
impurity spin: underscreening, with an asymptotically free effective impurity spin, and exact screening with a finite 
magnetic susceptibility of the impurity. Overscreened critical behavior is absent for this correlated electron model. As 
required by integrability, our model supports only forward scattering of the electrons off the impurity, and hence does 
not produce localization of the electrons. In contrast, more realistic ID magnetic impurity models must contain also 
back scattering terms which are expected to renormalize to an effective boundary potential, in addition to producing 
an impurity thermal response sensitive to the electron-electron interaction |20| . Yet, our model dramatically shows 
how correlations in the bulk make the magnetic impurity drive both spin and charge fluctuations. This property is 
expected to be generic and our exact solution provides a detailed picture of its possible realization. 
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